Effect of interfacial Cr on magnetoelectricity of Fe2/CrO2/BaTiO3(001) 
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On the basis of first-principles calculations we study the effect of interfacial Cr on the magne- 
toelectric properties of a composite multiferroic Fe_L/BaTiO3(001), with the Fe thickness L <2 
monolayers. The use of the Cr02-terminated interface instead of Ti02 niay significantly enhance 
magnetoelectricity in the system, showing an unexpected change in magnetization induced by the 
electric polarization reversal. In the case of L =2, for instance, the magnetic order of the Fe bilayer 
can be switched from nearly zero ferrimagnetic to ferromagnetic upon polarization reversal. 

PACS numbers: 31.15.Ae, 68.47.Gh, 73.20.At, 77.84.Lf, 77.80.Fm 
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I. INTRODUCTION 

The occurrence of ferroelectricity and ferromagnetism 
in the same phase of a so called multiferroic^ (MF) ma- 
terial allows both a switchable electric polarization, P, 
and a switchable magnetization M. More precisely, when 
an applied electric field displaces the magnetic ions of 
the multiferroic this affects the magnetic exchange cou- 
pling or, vice versa, the external magnetic field, H, in- 
duces Pi ^ aijHj, where a^ is the magnetoelectric 
(ME) tensor and («,j) = x,y,z. When a is sufficiently 
strong this phenomenon may allow to store information 
in nanometer-sized memories with four logic statea^Ti^. 

The classification of multiferroics is based on differ- 
ent mechanisms of induced polarity^ . The type-I class of 
multiferroics contains numerous perovskitelike materials 
in which P appears at higher temperatures than mag- 
netism. In these materials, P and M weakly interact 
with each other and, therefore, a is marginal there. In 
type-II MF, such as TbMnOa, ferroelectricity is driven by 
the electronic order degrees related to a spin-orbit mecha- 
nism in conjunction mostly with the spin-spiral magnetic 
arrangement via the Dzyaloshinskii-Moria antisymmetric 
exchange. The latter creates P ~ ry x [Si x Sj], where 
Tij is the vector connecting neighboring spins Si and S j . 
Some of the type-II MFs may disclose a relatively large 
ME coupling. However, their ferroelectricity is caused by 
a particular type of magnetic order, which exists only at 
low temperature and which is predominantly antiferro- 
magnetic. 

Studies based on density functional theory (DFT) have 
significantly contributed to this rapidly developing field 
of multiferroics^. For instance, calculations from first- 
principles predict that the ME effect appears when a meV 
voltage is applied across the interface between the two 
unlike terminations, such as SrRuOa/SrTiOy^. The in- 
terface ME effect might be intrinsically enhanced by the 
use of material with high spin polarization. Indeed, a 
more robust scenario of magnetoelectricity occurs in epi- 
taxially grown two-phase MF consisting of ferroelectric 
and ferromagnetic components. Ab initio calculations 
suggest that chemical bonding at the Fe/BaTiO3(001) 



interface is the source of strong ME coupling^' ^". More- 
over, for the two opposite directions of P (Pj, and Pf), 
there are rather noticeable differences of 0.1-0. 2/is in the 
magnetic moments of interfacial Fe and Ti. This is a very 
promising phenomenon, which is entirely confined to the 
ferroelectric/ferromagnetic interface. The interface ME 
effect^ defines the change in M at the coercive field Ec- 



fioAM « aEc 
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For Fe/BaTiO3(001), the estimatediS a of ~2xl0-i° G 
cm^/V is two orders of magnitude larger than that pre- 
dicted for SrRuOg/SrTiOa. 

Currently, ab initio calculations which explore the 
trends and basic physics of magnetoelectrics, go ahead 
of experiment. For a single Fe monolayer (ML) 
on BaTiO3(001), DFT predicts that perpendicular 
anisotropy is favored to in-plane anisotropy by 0.7 meV 
(0.5 meV) per Fe atom for P^ {P^)^. Although the spin 
reorientation transition under switching of P is not found 
from first principles, the ME coupling alters the mag- 
netocrystalline anisotropy energy by ~50%. The mag- 
netic order of Fe/BaTiOa can be tuned by the Fe layer 
thickness to almost zero-M ferrimagnetic upon deposi- 
tion of a second Fe MlJ^. Ferromagnetic order is re- 
stored for the Fe films thicker than 3 ML where the shape 
anisotropy energy favors in-plane alignment of M^^. Epi- 
taxial growth of the two-phase MF thin films of high 
quality continues to be very challenging. A 30-nm thick 
Fe(OOl) film has been grown recently on a ferroelectric 
BaTiO3(001) substrate^?.^ p^^. ^-^^^ composite MF, the 
trends of magnetic anisotropy are in good agreement 
with the corresponding ab initio calculations—'-^. Un- 
til recently, the DFT studies of the interface ME cou- 
pling were focused on chemically perfect films and su- 
perlattices with no impurities. Modeling the two dif- 
ferent Fe3O4/TiO2/BaTiO3(001) interfaces, within the 
DFT, Niranjan et alr^ have found that ME couphng 
is stronger for the 0-deficient type of the Fe304 inter- 
face. Therefore, the presence of extra oxygen or oxygen 
vacancies at the biferroic interface plays an important 
role. The effect of iron oxidation on the ME coupling of 
Fe/ATiO3(001) (^=Ba, Pb) was simulatedii from first 
principles for oxygen coverages ranged between 0.5 and 



2.0 adsorbed O atom per Fe atom. The calculations sug- 
gest that the magnetic properties of the Fe monolayer are 
gradually degraded with increasing O coverage. However, 
the change in magnetization which is induced by the P 
reversal remains robust. Thus, the surface oxidation of 
composite MFs cannot destroy their potentially switch- 
able magnetoelectricity. 

It is well known that both the magnetic order of Fe- 
films and the related magnetic anisotropy are very sen- 
sitive to the presence of some other 3d elements. The 
alloying effect may result in important changes in mag- 
netoelectricity and therefore, the DFT based modelling of 
chemical order in composite multiferroics would be use- 
ful. The effect of Fe-Co alloying on magnetoelectricity 
of thin-film Fe/BaTiO3(001) has been studied recentlyi^ 
from first principles using the coherent-potential approx- 
imation to DFT. It was found that the presence of >0.25 
Co at.% per Fe atom stabilizes the ferromagnetic order in 
the two-ML thick and magnetically soft Fe-films. In this 
work, we investigate the ME coupling in the 1-ML and 2- 
ML thick Fe on BaTiO3(001) (BTO), with a Cr02 inter- 
facial layer instead of Ti02. Chromium dioxide (Cr02) is 
an experimentally proven half metal, which shows a Curie 
temperature of 392 K and which possesses the largest 
spin polarization so far reported for this class of materi- 
als. As a consequence of the half-metallic feature of Cr02 
, the occupied Cr 3d bands are fully spin polarized, lead- 
ing to the spin moment of 2 fj,B per formula unit. Now 
we explore whether such a Cr02-terminated interface of 
BTO enhances the ME coupling in FeL/BTO(001). In 
particular, for L = 2 we observe a dramatic change of 
magnetization in the topmost Fe ML under polarization 
reversal. 



II. METHOD 

To model the FeL/BaTiO3(001) biferroic system 
within a slab geometry we used a 5-unit-cell (~2-nm) 
thick BTO supercell covered by an Fe monolayer or Fe bi- 
layer (L =1,2). A 2-nm- vacuum layer separates the slabs 
along [001]. For tetragonal BTO the equilibrium lattice 
parameters a =3.943 A and c/a =1.013 were used. The 
Fe positions and atomic positions of the two top BTO 
unit cells were relaxed. In ferroelectric BTO, the cations 
and O of the alternating BaO and Ti02 layers are dis- 
placed against each other in the [001] direction. This 
leads to spontaneous polarization along [001]. Here we 
model a dually polar ferroelectric. If the BTO cations 
are placed above O in the supercell then the negative in- 
tralayer displacements 6 — {zq — Zcation) < form the 
P state pointing parallel to the surface normal (P-j-) and, 
vice versa, the state P^ means that the S > 0. Before re- 
laxation, the S values of 0.082 A and 0.086 A were chosen 
in the Ti02 and BaO layers, respectively^^. The Ti02- 
terminated type of the BTO interface was energetically 
preferrediS. In this work, we substitute an interfacial Ti 
by Cr and added one or two ML of iron on the Cr02- 
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FIG. 1: Top layers of the (Fe2)L=2/CrO2/BaTiO3(001) unit 
cell are shown as side- and top view in the panels (a) and (b) , 
respectively. In panel (b), the interfacial Fe atoms above O 
are marked as "Feo" while the Fe atoms of the second ML, 
which sit above the perovskite cations are marked with the 
labels "Fesa" and "Fex". 



terminated BTO(OOl). The Fe adatoms of the first ML 
relax atop oxygen^'', while the Fe atoms of the second 
ML find their relaxed positions above the Ba and X=Cr 
sites. In Figure [U we plot the side- and top view of re- 
laxed FeL/CrO2/BTO(001) for the case of L = 2. The 
positions of Fe above O, Ba and X are indicated by the 
corresponding labels in the panel (b). 

In this DFT based study we used the Vienna Ab 
initio Simulation Package (VASP)i2r2S within the lo- 
cal spin-density approximation. The electron-ion in- 
teractions were described by projector-augmented wave 
(PAW) pseudopotentials^^, and the electronic wave func- 
tions were represented by plane waves with a cutoff en- 
ergy of 650 eV. For ionic relaxation the 8x8x4 k-point 
Monkhorst-Pack^^ mesh was used. The ionic relaxation 
was performed until the forces were less than 1 x 10^'^ 
eV/A. To calculate the electronic density of states (DOS) 
we used the 30 x 30 x 15 k-point mesh. For each com- 
pletely relaxed atomic configuration we performed the 
spin-polarized calculations starting form the ferromag- 
netic (FM) or, alternatively, from the antiferromagnetic 
(AFM) configuration in the Fe layers. The induced mag- 
netization of the XO2 interface was as well investigated. 



III. RESULTS AND DISCUSSION 

Much effort has been recently put to show that the 
electric field-induced reversal of P is able to vary the 
easy direction of magnetization in magnetically soft 
Coo.gFeo.i— and Nio.78Feo.22 permalloy^'^ attached to 
thin film of multiferroic BiFeOa (BFO) or, alternatively, 
to a single crystal of BiFeOs. There is a problem, how- 



ever, to form a ferroelectric single domain in the (001) 
plane of BFO. As a result, the magnetization of permal- 
loy could not be completely switched. We suggest that 
BTO is a more promising material for switching M by 
an electric field in the FM layer. In our study of the 1- 
ML-thick Fe-electrode material deposited on BTO(OOl) 
we find that the two systems: FeL=i/Ti02/BT0 and 
FeL=i/Cr02/BT0 are both ferromagnetically ordered, 
while the ME coupling coefficient increases from a = 
2.1 X IQ-io G cmVV in FeL=i/Ti02/BT0 to the value 
of 7.2 X 10^1° G cmVV at the CrOz interface. Eq.(l) was 
used to estimate a. In the case of Fe bilayer, the magnetic 
order changes dramatically. The FeL=2/Ti02/BTO sys- 
tem is almost zero-M ferrimagnetic for the both P states. 
Contrarily, FeL=2/Cr02/BTO changes its magnetic or- 
der from AFM to FM when the substrate polarization is 
switched from Pj^ to P-^, resulting in a ME coupling coef- 
ficient of a = 6 X 10~® G cm^/V. Below we concentrate 
mainly on the case of L = 2. 
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FIG. 2: Intralayer displacements 5 = zq — Zcauon 
(in A) calculated for several top perovskite layers of 
Fez,/XO2/BaTiO3(001) (L = 1, 2, X=Ti,Cr and P = Pj,, P^). 
The interfacial layer and layers beneath are denoted by I, I-l, 
1-2, etc. 



Structural relaxation 



Figure [H shows the perovskite intralayer displacement 
between oxygen and cations along [001], S = zq — Zcation, 
obtained after relaxation of FeL/XO2/BTO(001) (L = 
1,2, X=Ti,Cr and P = P;,Pt). The interfacial layer 
and layers beneath are denoted in Fig. [H by I, I-l, 1-2, 
etc. The asymmetry of 5 seen between Pj, and P-\ for the 
layers I, I-l and 1-2 as well as the magnitude of i5, which 
gradually decreases towards the interface, both mimic the 
effect of the depolarizing field and its screening. It should 
be noted that the state P^ is energetically preferred com- 
pared to Pf. For that reason the depolarization effect 
is rather strong for P^ as shown in Fig (21 For P^, the 
value of 5 is stable beneath the interface, namely, be- 
tween the layers I-l and 1-3 and, therefore, the reduction 
of i5 becomes crucial at the interface only. It turns out 
that interfacial Cr02 obeys marginal 5, which value de- 
creases when the second Fe ML is added. For P-j-, the 
effect of X=Cr on 5 is more pronounced. For instance, 
when L = 2 and P = P-^ the presence of Cr changes the 
sign of 6 in layer I. 

In Figure [31 we plot the relaxed distances between in- 
terfacial Fe and O atoms of XO2 (X = Ti, Cr). It has 
been previously found from first principles that the Ti02 
termination of BTO (001) is energetically preferred^^. 
When the first Fe ML is deposited on TiO2/BTO(001) 
the Fe atoms find their relaxed positions above O^- at 
the distance dj+ij «1.78 A as shown in the left panel 
of Fig. [21 Thus, Fe and O form a strong and relatively 
short chemical bond at the interface. Our calculations 
demonstrate that di+ij may increase by '^5 % when 
the second Fe ML is added. The polarization reversal 
shows no effect on dj+ij. For L = 1 and the Cr02- 
interface, we find that the corresponding di+ij «1.7 A 
is significantly reduced compared to the Fe/Ti02/BT0 
systems. When the Fe-(I-|-2) layer is added for X=Cr and 
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FIG. 3: Interlayer distances di+i./ between the Fe adlayer I-l-l 
and interfacial O are shown versus the Fe thickness L = 1, 2 
in the left panel. For L = 2, the distances d/+2,/+i between 
the topmost Fex/Fesa sites and the Fe-(I-l-l) layer are shown 
in the right panel. Each system is dually polar. 



P-]-, the separation between Fe and O is increased to the 
corresponding X=Ti value. For the opposite polarization 
Pj, and L = 2, a ~5 %-increase of dj^ij was obtained. 
The latter result suggests a very promissing scenario of 
magnetoelectricity in the FeL/Cr02/BT0 system with 
L = 2. Since the Fe-(I-|-2) atoms of the second layer are 
inevitably placed above the perovskite cations, the corre- 
sponding Fex and Fe^a sites are nonequivalent as shown 
in Fig. 1. In Fig. ^h) we plot the relaxed interlayer sep- 
aration di+2.i+i between the Fe layers I-I-2 and I-Fl for 
the case of L = 2. In general, the presence of Cr at the 
interface makes di+2,1+1 larger compared to the reference 
FeL/Ti02/BT0 system but, most importantly, di^2,i+i 
is not changed upon P reversal, exept for a 3 %-increase 
at the FcTi site. 



B. Electronic and magnetic properties 

Fig. m shows the site-projected DOS of paraelectric 
cubic BaTiOa together with the DOS of hypothetic cu- 
bic BaCrOa. The two perovskites were calculated using 
the same lattice parameter a = 3.943 A. For BTO we ob- 
tained an insulating band gap of ~2 eV, which is typically 
underestimated within the local density approximation. 
The conduction band of BTO is formed mainly by the 



in a = 7.2x10 ^° G cm^/V, which is three times larger 
than the ME effect of FeL=i/Ti02/BT0. 
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FIG. 4: The total and site-projected DOS of cubic BaTiOs 
and hypothetic BaCrOs calculated using the same lattice pa- 
rameter a — 3.943 A. 

Ti 3d states whereas the upper valence band is largerly 
composed by the O 2p states. In BaCrOs, the DOS is 
typically metallic while the 3d states of Cr dominate near 
the Fermi level, Ep ■ There is a marginal pseudogap seen 
at -1.5 eV below Ep. Therefore, one can expect relatively 
strong metallization at the Fe/Cr02 interface compared 
to Fe/TiOa. 

In FeL=i/Ti02/BT0 the FM order is energetically 
favorable against the AFM solution by 0.7 eV/cell 
(0.75 eV/cell) for P-j- (P^). Here, the Fe and O magnetic 
moments are aligned parallelly whereas the Ti magnetic 
moment, originating from hybridization of the Ti 3d and 
Fe 3d minority states^, is antiparallelly aligned. All mag- 
netic moments of the system are collected in Table HI The 
polarization reversal from P^ to P-]- yields the magnetiza- 
tion change |AM| = 0.028 /is/cell which formally results 
in the ME coupling of 2.1x10-1° G cmVV. When Cr 
substitutes Ti at the interface, the lowest-energy config- 
uration remains ferromagnetic. However, the negative 
magnetic moment of ~2 ^b, induced on Cr, is much 
larger than mxi- For interfacial oxygen the calculated 
magnetic moment is about 0.1 /is. This value as well 
as mcr are in a good agreement with the experimental 
data of bulk Cr02^. Due to the large and negative Cr 
magnetic moment, the total magnetization of the system 
FeL=i/Cr02/BT0 is reduced by ss 2 /zb in comparison to 
that of FeL=i/Ti02/BT0. Although mcr is moderately 
changed by P reversal the corresponding |AAf | results 




FIG. 5: Energy difference AE — Eafm — Efm between the 
AFM and FM configurations of FeL/X02/BT0 (X=Cr,Ti 
and L = 1,2) is normalized per Fe atom. 



The second Fe ML deposited on the Ti02-terminated 
BTO(OOl) interface causes a specific case. There are two 
inequivalent 1+2 sites situated atop Ba and Ti, respec- 
tively, which are labelled by Yesa and Yex in Fig. 1. 
The different magnetic moments refiect the neighbour- 
hood of these atoms such as their atomic volumes and 
hybridization of the electronic states. Let us consider, 
first, the case of X=Ti. The value of mpe in the layer 
I4-1 is almost quenched while the two sizable moments 
in the surface layer I-|-2 are antiparallelly aligned. This 
results in M ^ for FeL=2/TiO2/BTO(001). In the 
case of the Fe bilayer on the Cr02-terminated BTO, the 
lowest-energy configuration is antiferromagnetic for P^ 
and becomes ferromagnetic for P-|-. For this polarization 
the Fe magnetic moments in the layer I+l are far below 
their bulk value but the two Fe-(I-|-2) magnetic moments, 
which are ferromagnetically aligned to each other, con- 
tribute significantly to the total M . We estimate that 
the total magnetic moment of the system changes from 
M <0.3 ^B to >8 ^B per unit cell area upon polariza- 
tion reversal. Thus, the polarization reversal produces for 
X=Cr the effect of switchable magnetization. In Fig. [SJ 
the difference in energy, AE = Eafm — EpM , calculated 
between the AFM and FM configurations and normal- 
ized per Fe atom, is plotted. For X=Cr the 2-ML-thick 
Fe film represents a specific case of a magnetically soft 
system at fixed P. Nevertheless, any magnetic switch 
upon P reversal requires an energy which exceeds the 
coercive field value of BTO. 

To illustrate the interface ME coupling mechanism, 
we plot in Fig. [5] and Fig. [7] the spin density imbal- 
ance, (n+(r)— n^(r)), obtained under P-reversal near the 
interface of FeL=2/Ti02/BTO and FeL=2/Cr02/BTO, 
respectively. The (100) plane cutting through the X 
and O interfacial sites shows where the largest changes 
of the spin density occur and, hence, from where the 



TABLE I: Local magnetic moments (in /is) calculated for the two Fe adlayers labeled by I+l and 1+2 and interfacial X (X = 
Cr, Ti) and O of Fe_t/XO2/BaTiO3(001) {L = 1, 2). In the topmost Fe layer 1+2, there are two nonequivalent sites denoted as 
FeBa and Fex- The total magnetization Altot includes the contributions from the interstitials. The energy difference between 
the AFM and FM configurations calculated for each system at P =(-Pf i Pi) is shown in eV per cell. 



Site Layer 


(Fe2)L=] 


/CrOa/BTO 


(Fe2)L. 


=2/Cr02/BTO 


(Fe2)L=] 


/Ti02/BT0 


(Fea)!,. 


=2/Ti02/BTO 




Pt 


n 


Pt 


n 


P^ 


Pi 


P, 


n 


FeBa (1+2) 


— 


— 


+2.00 


+2.24 


— 


— 


+2.41 


+2.36 


Fex (1+2) 


— 


— 


+2.41 


-2.61 


— 


— 


-2.46 


-2.36 


Feo (I+l) 


+2.72 


+2.75 


+0.86 


+0.44 


+2.83 


+2.81 


-0.03 


+0.00 


X (I) 


-2.10 


-2.00 


+1.79 


-0.11 


-0.30 


-0.22 


0.00 


+0.01 


(I) 


+0.10 


+0.11 


-0.03 


0.00 


+0.09 


+0.08 


-0.01 


-0.01 


Eafm - EpMieY) 


+0.65 


+0.70 


+0.01 


-0.02 


+0.69 


+0.75 


-0.12 


-0.12 


Mtot (^s) 


+3.86 


+3.96 


+8.28 


+0.28 


+5.87 


+5.84 


+0.02 


-0.02 


a (G cmVV) 


7.20 


X 10-^" 


5.99 X 10"** 


2.08 


X 10-^" 


3.05x10-'" 



ME effect arises. Each of the four panels of Fig. |BH7] 
shows the local magnetization density calculated at fixed 
P ={Pi,P^). These are shown for the two possible mag- 
netic configurations which are either FM or AFM. For 
X=Ti, both the P^- and P|,-poled states are antiferro- 
magnetically ordered, as shown in the panels (b) and (d) 
of Fig. [SI The two results are similar to each other. The 
largest negatively charged areas are seen around Fex- 
(1+2) while the n+- charged areas around the second Fe 
site of this layer are not shown in Fig. [HI All other sites 
of FeL=2/Ti02/BTO including Fe-(I+1) indicate very 
small magnetic moments. Inspecting the spin density 
imbalance seen in FiglTJa) and FiglT^d) for the two ener- 
getically preferred but oppositely poled configurations of 
X=Cr, we find many differences in the magnetic struc- 
ture. The panel (a) shows the ferromagnetically ordered 
state P-j- where the Fe and Cr atoms form rather spa- 
cious regions of positive spin density n"*" while n" can 
be spotted around O, in the Fe interstitials and regions 
towards the surface. In the case of P^, the energetically 
favorable AFM configuration, shown in the panel (d), is 
similar to that of FeL=2/Ti02/BTO. Here, the large ar- 
eas around Fex-(I+2) and also around interfacial Cr are 
negatively charged. Besides, the p^-orbitals of interfacial 
O show their negative spin population resulting from hy- 
bridization with the 3d states of Fe-(I+1) whereas the 
O Pa; and Py orbitals, which form the bonds with the 
Cr 3d states, contribute to n"*". Regarding the Fe-(I+1) 
atoms of FeL=2/Cr02/BTO, Figl7i;d) shows that they 
contribute to n+ contrarily to the case of X=Ti. 

The site-projected and spin-resolved DOS calculated 
for FeL=2/X02/BTO are plotted in the two panels of 
Fig. [5] For each system, the solid (shaded) lines repre- 
sent the DOS curves in the Pf (P|_) state. The energet- 
ically preferable magnetic configurations are shown only 
in Fig. [H]for each direction of P. In general, the DOS of 
the interfacial XO2 layer is metallic for both systems. For 
L — 2 and X=Cr, however, the Cr Sd-DOS indicates rel- 
atively strong spin polarization at the Fermi level. This 
is not surprising since the DOS of hypothetical BaCrOa 
shows similar behavior, as shown in Fig. 4. When X=Ti, 
there is some insignificant presence of the Ti 3d states 
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FIG. 6: Spin density imbalance n"'"(r) — n (r) (in e/A ^) 
within the (100) plane cutting through the Ti atoms of the 
FeL=2/TiO2/BaTiO3(001) slab. The black vertical lines rep- 
resent the unit cell boundary. The two top (bottom) panels 
show the polarization Pi^ (Pi)- The panels (a) and (c) illus- 
trate the FM ordering obtained for Pi^ and Pj,, respectively, 
while the energetically preferable and nearly AFM configura- 
tions are shown in the panels (b) and (d). 



in the BTO band gap below Ep, which entirely results 
from the hybridization with the Fe 3d states of the layer 
I+l. Another major difference in the DOS seen in Fig. 
|H] for L = 2 comes from the magnetic ordering of Fex ■ 
For X=Ti the two Fe atoms in the topmost layer 1+2 
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FIG. 7: Spin density imbalance (in e/A ^) within 
the (100) plane cutting through the Cr atoms of 
Fei=2/CrO2/BaTiO3(001). The two top (bottom) panels 
show the Pf (Pi) states while the left (right) panels illustrate 
the FM (AFM) ordering. For P^ (Pi), the lowest energy so- 
lution is the FM (AFM) configuration shown in 'a' ('d'). 



are coupled antiferrimagnetically while the corresponding 
DOS curves show minor changes upon P reversal. When 
X=Cr the polarization reversal from P^ to the state P-t- 
supports (i) the ferromagnetic order in the layer 1+2, (ii) 
the relatively large magnetic moment TO(Feo)~0.9 fiB in 
the layer I+l and (iii) the r^2-fiB change of rricr which 
is aligned parallelly to the Fe magnetic moments. 

In Fig. [SI we plot the relative (in %) and absolute con- 
tributions (in hb) to AM = M{P^)—M{P^) coming from 
each magnetic species of FeL=2/X02/BTO. For the two 
biferroic interfaces studied here, the largest P-induced 
change of M comes from the Fe-(I-|-2) atoms. For X=Cr, 
however, the absolute value of AM approaches ^7 fj,B per 
unit cell. As result, the corresponding ME coupling coef- 
ficient increases significantly compared to that of X=Ti. 
We demonstrate that the case of L = 2 and X=Cr sta- 
bilizes the FM ordering in the system with P pointing 
upwards. Surprisingly, this is completely due to rather 
modest 5 % decrease oidj+ij under P reversal, as shown 
in Fig. 3. With decreasing dj+ij above Cr02, the FM 
order is developing in the system. More precisely, when 
the Fe-(I-l-l) magnetic moment becomes larger ferromag- 
netism is stabilized in layer (1+2). In the case of X=Ti, 
the interlayer separations (i/+2.7+i and dj+i./ are almost 
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FIG. 8: Spin-polarized and site-projected DOS cal- 
culated for the metal 3d- and O 2p-states near the 
inerface of (Fe2)L=2/CrO2/BaTiO3(001) (left) and 
(Fe2)L=2/TiO2/BaTiO3(001) (right). The two upper 

panels show the Fe 3d DOS of the surface layer 1+2 while the 
Fe-I+1 DOS are shown in the middle panels. The 3d-DOS of 
inerfacial cations Cr/Ti and oxygen 2p-DOS are plotted in 
the two lower panels. Solid lines and shaded areas represent 
the DOS for P^ and P^, respectively. 



the same upon the P reversal that prevents any crucial 
spin reorientation in the topmost Fe ML. 



IV. CONCLUSIONS 

In summary, we present an ab initio study of the 
effect of interfacial Cr on the strength of magneto- 
electric coupling seen at the interface of multiferroic 
FeL/CrO2/BaTiO3(001), with the Fe thickness L <2 
monolayers. We predict that a Cr02-terminated inter- 
face instead of Ti02 may significantly enhance magneto- 
electricity in the system. The most attractive scenario 
is, however, obtainded for the Fe bilayer where the mag- 
netic order changes from nearly zero-M ferrimagnetic to 
ferromagnetic upon polarization reversal in ferroelectric 
BaTiO3(001). 
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